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Abstract 



a" 

We consider the prospects for dark matter /energy unification in k-essence type 
(3J[)! cosmologies. General mappings are established between the k-essence scalar field, the 

hydrodynamic and braneworld descriptions. We develop an extension of the general 
relativistic dust model that incorporates the effects of both pressure and the associated 
acoustic horizon. Applying this to a tachyon model, we show that this inhomogeneous 
"variable Chaplygin gas" does evolve into a mixed system containing cold dark matter 
like gravitational condensate in significant quantities. Our methods can be applied to 
any dark energy model, as well as to mixtures of dark energy and traditional dark 
< 3\ ' matter. 
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1 Introduction 



The discovery of the accelerated Hubble expansion in the SNIa data combined with 
observations of the cosmic microwave background (CMB) [2J, [3] , has forced a profound shift 
in our cosmological paradigm. If one makes the conservative assumptions of the validity of 
Einstein's general relativity and the cosmological principle, one concludes that the universe 
is presently dominated by a component that violates the strong energy condition, dubbed 
dark energy (DE). Moreover, primordial nucleosynthesis constrains the fraction of closure 
density in baryons, Qb, to a few percent, while galactic rotation curves and cluster dynamics 
imply the existence of a nonbaryonic dark matter (DM) component with Qdm ^ (for a 
review see [1]). Currently, the best fit values are Qb = 0.04, f2 d^j = 0.22 and VLde = 0.74 
[3]. Thus it may be said that we have a firm theoretical understanding of only 4% of our 
universe. 

Pragmatically, the data can be accommodated by combining baryons with conventional 
cold dark matter (CDM) candidates and a simple cosmological constant A providing the 
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DE. This ACDM model, however, begs the question of why A is non-zero, but such that DM 
and DE are comparable today. The coincidence problem of the ACDM model is somewhat 
ameliorated in quintessence models which replace A by an evolving scalar field. However, like 
its predecessor, a quintessence-CDM model assumes that DM and DE are distinct entities. 
For a recent review of the most popular DM and DE models, see [5]. 

Another interpretation of this data is that DM/DE are different manifestations of a 
common structure. Speculations of this sort were initially made by Hu [6]. The first definite 
model of this type was proposed a few years ago [7J [9], based upon the Chaplygin gas, a 
perfect fluid obeying the equation of state 

P =--, ID 

which has been extensively studied for its mathematical properties [10J. The general class 
of models, in which a unification of DM and DE is achieved through a single entity, is often 
referred to as quartessence [III 112] . Among other scenarios of unification that have recently 
been suggested, interesting attempts are based on the so-called k-essence [TBI EH], a scalar 
field with noncanonical kinetic terms which was first introduced as a model for inflation |15j . 

The cosmological potential of equation ([1]) was first noted by Kamenshchik et al [TJ, who 
observed that integrating the energy conservation equation in a homogeneous model leads 
to 




p(a) = \A+-, (2) 



where a is the scale factor normalized to unity today and B an integration constant. Thus, 
the Chaplygin gas interpolates between matter, p ~ \^Ba~ 3 , p ~ 0, at high redshift and a 
cosmological constant like p ~ y/A ~ — p as a tends to infinity. The essence of the idea in [HIE] 
is simply that in an inhomogeneous universe, highly overdense regions (galaxies, clusters) 
have \w\ = \p/p\ <C 1 providing DM, whereas in underdense regions (voids) evolution drives 
p to its limiting value yA giving DE. 

Of particular interest is that the Chaplygin gas has an equivalent scalar field formulation 
[HI El [TO] • Considering the Lagrangian 

C = -y/Ay/l-X, (3) 

where 

X = , (4) 

equation (TjQ) is obtained by evaluating the stress-energy tensor T M „, and introducing = 
(p t/J/ / \I~X for the four- velocity and p = \fAj \J\ — X for the energy density. One recognizes C 
as a Lagrangian of the Born-Infeld type, familiar in the D-brane constructions of string/M 
theory [TB]. Geometrically, C describes space-time as the world- volume of a 3+1 brane in a 
4+1 bulk via the embedding coordinate X 5 [17] . 

To be able to claim that a field theoretical model actually achieves unification, one must 
be assured that initial perturbations can evolve into a deeply nonlinear regime to form a 
gravitational condensate of superparticles that can play the role of CDM. In [HI [9] this was 
inferred on the basis of the Zel'dovich approximation [T5]. In fact, for this issue, the usual 
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Zel'dovich approximation has the shortcoming that the effects of finite sound speed are 
neglected. 

All models that unify DM and DE face the problem of nonvanishing sound speed and the 
well-known Jeans instability. A fluid with a nonzero sound speed has a characteristic scale 
below which the pressure effectively opposes gravity. Hence the perturbations of the scale 
smaller than the sonic horizon will be prevented from growing. Soon after the appearance 
of [7J and [8], it was pointed out that the perturbative Chaplygin gas (for early work see 
|19j . and more recently [20]) is incompatible with the observed mass power spectrum [21] 
and microwave background [22]. Essentially these results follow from the adiabatic speed of 
sound 

T2 ( 5 ) 



2 _ 
° S dp 



which leads to a comoving acoustic horizon 



d s = fdt— . (6) 
J a 

The perturbations whose comoving size R is larger than d s grow as 5 = (p — p)/p ~ a- Once 
the perturbations enter the acoustic horizon, i.e., as soon as R < d s , they undergo damped 
oscillations. In the case of the Chaplygin gas we have d s ~ a 7 / 2 /H , where H is the present 
day value of the Hubble parameter, reaching Mpc scales already at redshifts of order 10. 
However, to reiterate a point made in [8], small perturbations alone are not the issue, since 
large density contrasts are required on galactic and cluster scales. As soon as 5 ~ 1 the 
linear perturbation theory cannot be trusted. An essentially nonperturbative approach is 
needed in order to investigate whether a significant fraction of initial density perturbations 
collapses in gravitationally bound structure - the condensate. If that happens the system 
evolves into a two-phase structure - a mixture of CDM in the form of condensate and DE in 
the form of uncondensed gas. 

The case, where the Chaplygin gas is mixed with CDM, has been considered in a number 
of papers [231 [2H 1251 [26l [271 12HI I2SI [30]. Here, the Chaplygin gas simply plays the role of 
DE. In keeping with the quartessence philosophy, it would be preferred if CDM could be 
replaced by droplets of Chaplygin gas condensate, as in [31 J . Homogeneous world models, 
containing a mixture of CDM and Chaplygin gas, have been successfully confronted with 
lensing statistics [231 121] as well as with supernova and other tests [251 I.26J. 

Another model, the so called "generalized Chaplygin gas" [32], has gained a wide popu- 
larity. The generalized Chaplygin gas is defined as [H [TJ [32] p = —A/p a with < a < 1 for 
stability and causality. As in the Chaplygin gas case, this equation of state has an equivalent 
field theory representation, the "generalized Born-Infeld theory" [321 133] . However, the as- 
sociated Lagrangian has no equivalent brane interpretation. The additional parameter does 
afford greater flexibility: e.g. for small a the sound horizon is d s ~ y/aa 2 /H Q , and thus by 
fine tuning a < 10 -5 , the data can be perturbatively accommodated [21] . Bean and Dore 
[27] and similarly Amendola et Al [28] have examined a mixture of CDM and the generalized 
Chaplygin gas against supernova, large-scale structure, and CMB constraints. They have 
demonstrated that a thorough likelihood analysis favors the limit a — > 0, i.e. the equivalent 
to the ACDM model. Both papers conclude that the standard Chaplygin gas is ruled out 
as a candidate for DE. However, analysis [301 E3] of the supernova data seems to indicate 
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that the generalized Chaplygin gas with a > 1 is favored over the a — > model and similar 
conclusions were drawn in [20]. But one should bear in mind that the generalized Chaplygin 
gas with a > 1 has a superluminal sound speed that violates causality [34]. For a different 
view on this issue see (35j [37] (see discussion in section [2]). 



The structure formation question, in respect of the Chaplygin gas, was decided in [38] . 

In fact, in the Newtonian approximation, we derived an extension of the spherical model 
[39] that incorporates nonlinearities in the density contrast 5, as well as the effects of the 
adiabatic speed of sound. Both are crucial, since for an overdensity we have c s < c s , where 
c s is the speed of sound of the background. Although small initial overdensities follow 
the expected perturbative evolution, for an initial <5i?(a; n ) exceeding a scale independent 
critical S c , it was shown that <5r(o) tends to infinity at finite redshift, signaling the formation 
of a bound structure or condensate. Unfortunately, it was further found that, when the 
required 5 C is folded with the spectrum of the initial density perturbations to obtain the 
collapse fraction, less than 1% of the Chaplygin gas ends up as condensate. Thus the simple 
Chaplygin gas is not viable due to frustrated structure formation. In effect, the model is 
a victim of the radiation dominated phase, which turns the Harrison- Zel'dovich spectrum 
6 k ~ k 1 / 2 to 6 k ~ k~ 3/2 at Rqeq ^ 26 Mpc. In a pure Chaplygin-gas universe without 
radiation there would inevitably be sufficient small scale power to drive condensation. 

One way to deal with the structure formation problem, is to assume entropy perturbations 
P, HO] such that the effective speed of sound vanishes^. In that picture we have 5p = 
c 2 5p — SA/p = even if c s ^ 0. But as we detail below, in a single field model it is precisely 
the adiabatic speed of sound that governs the evolution. Hence, entropy perturbations require 
the introduction of a second field on which A depends. Aside from negating the simplicity of 
the one-field model, some attempts at realizing the nonadiabatic scenario [HJ H2J H3] have 
convinced us that even if Sp = is arranged as an initial condition, it is all but impossible 
to maintain this condition in a realistic model for evolution. 

The failure of the simple Chaplygin gas does not exhaust all the possibilities for quart- 
essence. The Born-Infeld Lagrangian ([3]) is a special case of the string-theory inspired tachyon 
Lagrangian [HI H5] in which the constant \fA is replaced by a potential V((p) 



In turn, tachyon models are a particular case of k-essence [15] . The possibility of obtaining 
both DM and DE from the tachyon with inverse square potential has been speculated in [46] . 
More recently, it was noted [17] that, in a Friedmann- Robertson- Walker (FRW) model, the 
tachyon model is described by the equation of state (CD) in which the constant A is replaced by 
a function of the cosmological scale factor a, so the model was dubbed "variable Chaplygin 
gas". Related models have been examined in [18], H9], however, those either produce a larger 
d s than the simple Chaplygin gas [15] , or else need fine-tuning [49] J§ 

In this paper we develop a version of the spherical model for studying the evolution of 
density perturbations even into the fully nonlinear regime. Although similar in spirit to 

1 Note this "silent quartessence" is not different from [8J , where we tacitly neglected the effects of nonva- 
nishing c s ■ 

2 The tachyon model [35] gives d s ~ a 2 /H n . The two-potential model [53] yields d s ~ \fl — ha 2 /Ho, so it 
requires 1 — h < 1CP 5 like the generalized Chaplygin gas. Expanding in 1 — h, the second potential reveals 
itself to be dominantly a cosmological constant. 





(7) 
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[38] , the formalism here is completely relativistic, rather than Newtonian, and applicable to 
any k-essence model instead of being restricted to the simple Chaplygin gas. The one key 
element we carry over from [38] is an approximate method for treating the effects of pressure 
gradients - which is to say the adiabatic speed of sound - on the evolution. Our method 
is flexible enough and can be extended to deal with a mixture of DE and DM. A spherical 
model has been applied to DE/DM mixtures in [50J, however there the effects of pressure 
gradients were omitted. 

We apply our method to the preliminary analysis of a unifying model based on the 
tachyon type Lagrangian ((Zj) with a potential of the form 

V(<p) = V n tp 2n , (8) 

where n is a positive integer. In the regime where structure function takes place, we show 
that this model effectively behaves as the variable Chaplygin gas with A(a) ~ a 6n with 
n = 1(2) for a quadratic (quartic) potential. As a result, the much smaller acoustic horizon 
d s ~ a ( - 7//2+3ra - ) / Hq enhances condensate formation by two orders of magnitude over the simple 
Chaplygin gas (n = 0). Hence this type of model may salvage the quartessence scenario. 

The remainder of this paper is organized as follows. In section [2] we reformulate k-essence 
type models in a way that allows us to deal with large density inhomogeneities. In section 
3 we develop the spherical model approximation that closes the system of equations. These 
two sections are completely general and stand alone. Numerical results, in section 4, are 
presented for positive power-law potentials and contrasted with the simple Chaplygin gas. 
Our conclusions and outlook are given in section [5j Finally, in |AJ we derive the adiabatic 
speed of sound for a general k-essence fluid, and in [Bj we give a brief description of the 
tachyon model from the braneworld perspective. 



2 K-essentials 

A minimally coupled k-essence model [13 [51], is described by 

R 



S 



d x 



-g 



16nG 



+ C(<p,X) 



(9) 



where C is the most general Lagrangian, which depends on a single scalar field ip of dimension 
m _1 , and on the dimensionless quantity X defined in For X > that holds in a 

cosmological setting, the energy momentum tensor obtained from takes the perfect fluid 
form, 

= 2£ x WW _ ^fiV = (P + P)u^u v -pg^v, (10) 
with Cx denoting dC/dX, and 4-velocity 



u, 



7} 



X 



(11) 



where rj is +1 or —1 according to whether tp t Q is positive or negative, i.e., the sign of w M is 
chosen so uq is positive. The associated hydrodynamic quantities are 



p = 



(12) 
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p = 2XC x (<p,X)-£(<p,X). (13) 
Two general conditions can be placed upon the functional dependence of £, 

£x > (14) 

and 

Cxx > 0. (15) 

The first condition (|T4"1) stems from the null energy condition, T^n^n" > for all light like 
vectors required for stability [52]. For a perfect fluid we have T liv n^n u = (p + p)(ii /t n /i ) 2 , 
thus p + p = 2XCx > 0, and owing to X > 0, we arrive at ffT4"|) . 

The second condition f|T5|) arises from restrictions on the speed of sound. Observing that 
f|T3|) allows us to view X as a function of p and the adiabatic speed of sound 



2 = 

° s ~ dp 



dp 



s/n 



(16) 

as shown in Appendix 0, coincides with the so called effective speed of sound 



s Px £x + 2X£ X x 

obtained in a different way in [51]. For hydrodynamic stability we require c~ > 0. In 
addition, it seems physically reasonable to require c\ < 1 in order to avoid possible problems 
with causality [31]. Causality violation in relation to superluminal sound propagation has 
been the subject of a recent debate [3H ESI ESI S3 E3]- It has been shown [351 E3] that if 
k-essence is to solve the coincidence problem there must be an epoch when perturbations 
in the k-essence field propagate faster than light. Hence, on the basis of causality it has 
been argued [53] that k-essence models which solve the coincidence problem are ruled out as 
physically realistic candidates for DE. In contrast, it has been demonstrated [35], [36J, [37] that 
superluminal sound speed propagation in generic k-essence models does not necessarily lead 
to causality violation and hence, the k-essence theories may still be legitimate candidates for 
DE. As the coincidence problem is not the issue in DE/DM unification models we stick to 
< c 2 s < 1. In view of X > and flU]) equation (JT5J) follows. 
Formally, one may proceed by solving the ip field equation 

(2£ x <ryv) ;At - = o, (18) 

with £ v denoting d£/d<p, in conjunction with Einstein's equations to obtain p and p. How- 
ever, it proves more useful to pursue the hydrodynamic picture. For a perfect fluid the 
conservation equation 

T^. v = (19) 
yields, as its longitudinal part u^T^ v - u = 0, the continuity equation 

p + 3H(p + p) =0, (20) 
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and, as its transverse part, the Euler equation 



where we define 
The tensor 



= , (21) 

P + P 



m = u v , v ; p = u v p, v ; = u v vP. v . (22) 



hfu, = g^u - u^u u (23) 

is a projector onto the three-space orthogonal to u M . The quantity Tt is the local Hubble 
parameter. Overdots indicate the proper time derivative. 

Now, using (fTTl) - (|T3l) the <p field equation ([TBI) can be expressed as 

(v^A -^ = 0, (24) 



or 



p + 3H(p + p) + (<p-r}Vx)£ v = 0. (25) 

where, as before, r] = sgn(<p )- Hence, for £ ip = (purely kinetic k-essence), the <p field 
equation (TT51) is equivalent to f[2"Uj) . In the general case, i.e. for £ ip ^ 0, equation ffTSl) 
together with (TSUI) implies 

2 = ^,P), (26) 

provided ffl3|) is invertible. 

Next, we observe that Euler 's equation can be written in various forms. Equation 

*„ = ^ (27) 

follows directly from (|TT|) - (fT3"|) and (|22p . With X a function of <p and p, the pressure p also 
becomes a function of <p and p. Thus by ( TTTj) . ffl6l) and fTSBl we find 



p + p 

This is a simple demonstration of the observation made earlier, that in a single component 
system it is the adiabatic (rather than effective) speed of sound that controls evolutional 
Rather than specifying £ directly, one may choose p = p(ip,p) and then find <? s from (TT2"j) . 
Up to an overall multiplicative integration function of ip only (which in turn can be absorbed 
in a reparameterization of (p itself) equations (j2"7|) and ff28l) imply 

X{<p t p) = Bxp(2f f^j , (29) 

which may be used in (1261) as an evolution equation for ip, while (120]) is an evolution equation 
for p. Further, equation (|29|) can be formally inverted to give p as a function of ip and X, 
thus allowing the construction of £ = p (ip, p(ip, X)). An example of this will be given in 
section 4. However, first we need to close the system of evolution equations. 



3 It seems to us that this point is often confused in the literature. 
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3 The Spherical Model 



Since the 4-velocity f TTTT) is derived from a potential, the associated rotation tensor vanishes 
identically. The Raychaudhuri equation for the velocity congruence assumes a simple form 

m + m 2 + a^ v + uV = u% . (30) 

with the shear tensor defined as 

<V = h a ^h 13 v v,( a .p) - Hhfa, . (31) 

We thus obtain an evolution equation for Ji that appears in (120]) . sourced by gravity through 
the Ricci tensor IZ^ and by both shear and the divergence of the acceleration If = 0, 
as for dust, equations ( l20l) and ( l30l) . together with Einstein's equations for R^ u comprise 
the spherical model [39]. However, we are not interested in dust, since generally a' 1 ^ as 
given by Euler's equation ( 1281) . Indeed, this term is responsible for the Jeans phenomenon in 
perturbation theory. One is only allowed to neglect w M in the long wavelength limit, where 
everything clusters, but one has no realistic information about the small (i.e. subhorizon) 
scales. 

The spherical top-hat profile is often invoked to justify neglecting the acceleration term 
(see e.g. ]EU\ and references therein). In fact, this leads to infinite pressure forces on the 
bubble boundary. Even if suitably regularized, the influence of these large forces on the 
bubble evolution is never accounted for. Needless to say, this makes the reliability of the 
inferences highly problematic, unless one invokes entropy perturbations again. 

For a one-component model, the Raychaudhury equation ( 1301) combines with Einstein's 
equations to 

m + m 2 + a 2 + 4vrG(p + 3p) = ( ^-fA . ( 32 ) 

V P + P )■„ 

with 

a 2 = a^. (33) 
In general, the 4-velocity can be decomposed as [51] 



u » = (if" + v") /Vl - v 2 , (34) 

where U 11 = 5q / ^/goo is the 4-velocity of fiducial observers at rest in the coordinate system, 
and is spacelike, with v^v^ = —v 2 and £/ M t> M = 0. In a mult i- component model (e.g. dark 
energy plus CDM, or a mixture of condensed and uncondensed k-essence) one can, in the first 
approximation, neglect the relative peculiar velocities. Then there will be a pair of equations 
fl2"U|) and (13"2"|) for each component with p + 3p replaced by a sum over all components and 
each component may be treated in comoving coordinates. In comoving coordinates v M = 
and 

= y/gooS°; = —=5q\ h 00 = h i = 0; hy = </y. (35) 
Then the nonvanishing components of the shear tensor are o"^- and from fl3Tl) it follows 

Cij = - Tihij , (36) 
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Assuming 

ds 2 = N 2 dt 2 — •jijdxidxj (37) 
we find ^ 

°ij = ~^ dt ^i + ■ ( 38 ) 
In spherically symmetric spacetime it is convenient to write the metric in the form 

ds 2 = N(t, rfdt 2 - b(t, r) 2 (dr 2 + r 2 f(t, r)dtt 2 ) (39) 

where N(t,r) is the lapse function, b(t,r) is the local expansion scale, and f(t,r) describes 
the departure from the flat space for which f — 1. We assume that N, a, and / are arbitrary 
functions of t and r which are regular and different from zero at r = 0. Then, the local 
Hubble paprameter and the shear are given by 

and 

In addition to the spherical symmetry we also require an FRW spatially flat asymptotic 
geometry, i.e., for r — » oo we demand 

JV-1; /->1; 6->a(t). (42) 

Here a denotes the usual expansion scale. 

The righthand side of (1321) is difficult to treat in full generality. As in [38], we apply the 
"local approximation" to it: The density contrast 5 = (p — p)/p is assumed to be of fixed 
Gaussian shape with comoving size R, but time-dependent amplitude, so that 

p(t,r)=p(t)[l+5 R (t)e- r2 ^ 2R \ (43) 

and the spatial derivatives are evaluated at the origin. This is in keeping with the spirit of 
the spherical model, where each region is treated as independent. 

Since dip = at r = 0, naturally diN = and d{b = at r = 0. Hence, 

N(t,r) = N(t,0)(l + O(r 2 )); b(t, r) = b(t, 0)(1 + 0(r 2 )). (44) 

Besides, one finds j\o — ► as r — > which follows from Einstein's equation G 1 o = 

*M$-%)M* + $) + Hl-W)-°- (45) 

Since the first three terms on the lefthand side vanish as r — > the last term can vanish if 
and only if f t o/f = at r = 0. By ( 14T1) the shear scalar a vanishes at the origin. 
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From now on we denote by Tl, 6, and N the correspobding functions of t and r evaluated 
at r = 0, i.e., H = Hit, 0), b = b(t, 0) and N = N(t, 0). According to (00]), the local Hubble 
parameter at the origin is related to the local expansion scale as 

Evaluating (1521) at the origin yields our working approximation to the Raychaudhuri equa- 
tion, i.e. we obtain 

Nlt +n + — {p + 3p) -b*Ri( P + p)- (47) 

The recommendations of this equation are that it extends the spherical dust model, by incor- 
porating both pressure and, via the speed of sound, the Jeans' phenomenon. In particular, 
it reproduces the linear theory with the identification k = y3/R for the wavenumber. 



4 Cosmological Tachyon Condensation 



We will now apply our formalism to a particular subclass of k-essence unification models 
described by (JTj). However, first it is useful to see how such models can be reconstructed 
using the methods of section 2. 

Violating the strong energy condition with positive p requires p < 0, while stability 
demands (? s = dp /dp > 0. These criteria are met by0 



P 



A{<p) 



A{ip) > 0, 



(4f 



for which 



* n a+l — 



a > 0. 



(49) 



Note that when the null energy condition is saturated, we have p a+l = A(ip), and that 
causality restricts a to a < 1. Using (|29|) . we arrive at 



X(<p,p) 



and the Lagrangian density of the scalar field 



2a/(l+a) 



c 



a/(l+a) 



X (l+a)/2a 



(50) 



(51) 



Only for a = 1 does one have (? s = 1 at the point where the null energy condition is saturated. 
Moreover, only for a = 1 can one obtain the tachyon Lagrange density 



c= -JA{<p)VT=x 



(52) 



J We do not consider the trivial generalization of adding a function of if alone to p. 
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which coincides with (jTj), identifying A(ip) = V(f) 2 . The equation of state is then given by 

P = -^, (53) 

P 

and the quantity X may be expressed as 



\2 



X(p^) = l-^^ = l-cl = l + w . (54) 
P 2 

Finally, only for a — 1, can the tachyon model be reinterpreted as a 3 + 1 brane, moving in 
a warped 4 + 1 spacetime, as shown in Appendix B. 

Equations , f[2"Ul) . f[2"B"j) . (JH)]) . and fT4T|) determine the evolution of the density contrast. 
However, as this set of equation is not complete, it must be supplemented by a similar set 
of equations for the background quantities p and H 

^ + 3(p + p) = 0, (55) 

— + H 2 + —(p + 3p) =0, (56) 

where p = p(p, (p). Due to (11) the field p in comoving coordinates is a function of time only, 
its gradient is always timelike or null, i.e., X > 0, and the perfect fluid description remains 
valid even in a deep non-linear regime. In comoving coordinates equation (26) reads 

(I)' = (57) 
In particular, in the asymptotic region r — > oo we have 

= (58) 

Equating fl57|) with fl58|) we find an expression for the local lapse function N in terms of X 



N = JX(<p,p)/X(<p,p). (59) 



Hence, the complete set of equations for p, H, (p, b, p, and H, consists of fl55l) . fl56l) . fl58|) . 
and 

% + WH(p + p) = 0, (61) 



~5~ 



(62) 
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where N is given by (|59|) . In this way, we have a system of six coupled ordinary differential 
equations that describes the evolution of both the background and the spherical inhomo- 
geneity. The definition of the Hubble parameter 



H 



1 da 
a dt 



(63) 



is used to express the evolution in terms of the background scale factor a. 

Here we restrict our attention to the power-law potential (JED- In the high density regime, 
where c 2 is small, we have X ~ 1, and (1261) can be integrated yielding ip ~ 2/(3H), where 
H ~ Ho^/Iia~ 3 / 2 , being the equivalent matter content at high redshift. Hence, as promised, 



Aitp) = V{<p) 



2 ~ a 6 ", 



which leads to a suppression of 10 at z = 9 for n = 1. 



To proceed we require a value for the constant V n in the potential (jSJ). Changing 14 
effects not only the pressure today but also the speed of sound and thereby the amount of 
structure formation. One must also be mindful that if there is a large amount of nonlinear 
structure formation then the single fluid description will break down at low redshift. As 
the main purpose of this paper is to investigate the evolution of inhomogeneities we will not 
pursue the exact fitting of the background evolution. Hence, rather than attempting to fit 
V n using the naive background model, we estimate V n as follows. We integrate (1261) with 



X = l + w(a) ~ 1 



as in a ACDM universe |48| and obtain 



(64) 



ip{a) 



arctan 




We then fix the pressure given by (J7J) to equal that of A at a = 1, i.e. 
yielding 



p n A = -v n ip(i) 2n 

3a n 



Vn 



8nG 



H, 



2(n+l) 




where 



a, 



arctan 




-2n 



(65) 

(66) 
(67) 

(68) 



so ao — 0.854, a.\ ~ 1.34, and — 2.09. With these values the naive background in our 
model reproduces the standard cosmology from decoupling up to the scales of about a = 0.8 
and fits the cosmology today only approximately. 

Using ( 1631) . we solve our differential equations with a starting from the initial dd ec = 
l/(#dec + 1) & t decoupling redshift z^ec = 1089 for a particular comoving size R. The initial 
conditions for the background are given by 



Po- 



tt 

2 dec 



Hir, = Hft 



n 

~~3 ' 

3 dec 



Pi 



3H\ r 



(69) 
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Figure 1: Evolution of 5 R (a) in the tachyon spherical model from a de c — 1/1090 for n — 2, 
R = 10 kpc, 5 R (a dcc ) =0.004 (solid) and 5 R (a dec ) =0.0055 (dashed). 



and for the initial inhomogeneity we take 

Pin = Pin(l + S- m ) , H m = H m (l - , (f in = <An = (70) 

where Q = 0.27 represents the effective dark matter fraction and 5i n = 8n(adec) is a variable 
initial density contrast, chosen arbitrarily for a particular R. 

Here it is worthwhile mentioning that the set of equations (155]) . (j581) . and (16D1) - 

preserves the condition X > and hence, the perfect fluid description remains valid 
throughout the nonlinear evolution described above. The evolution starts from a small 
initial value of <^i n and large p m so that initially X ~ 1 and we take initial <fi to be positive. 
According to (154|) . as </j increases and p decreases, X decreases up to a point where it becomes 
0. The first such point may be roughly at a between 0.1 and 1. At that point the sign of (p 
flips and X remains positive up to the next point at which X = 0. The sign of Lp flips again 
keeping X positive, and so on. Basically, this process continues ad infinitum never violating 
X > 0. Hence, the perfect fluid assumption is legitimate even in the deep non-linear regime. 

In figure [1] the representative case of evolution of two initial perturbations starting from 
decoupling for R = 10 kpc is shown for n = 2. The plots represent two distinct regimes: the 
growing mode or condensation (blue dashed line) and the damped oscillations (red solid line). 
In contrast to the linear theory, where for any R the acoustic horizon will eventually stop 5r 
from growing, irrespective of the initial value of the perturbation, here we have for an initial 
^(fldec) above a certain threshold S C (R), 5r(o) — > oo at finite a, just as in the dust model. 
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Figure 2: Initial value 5.R(adec) versus R for fl = 0.27 and h = 0.71. The threshold 6 C (R) is 
shown by the line separating the condensation regime from the damped oscillations regime. 
The solid line gives cr(R) calculated using the concordance model. 



Thus perturbations with Sn(adec) > S C {R) evolve into a nonlinear gravitational condensate 
that at low z behaves as pressureless super-particles. Conversely, for a sufficiently small 
<5R( a dec), the acoustic horizon can stop 5r(o) from growing; at low redshift the perturbations 
behave as expected from linear theory. Figure [2] shows how the threshold 5 C (R) divides the 
two regimes depending on the comoving scale R. 

The crucial question now is what fraction of the tachyon gas goes into condensate. In 
[31] it was shown that if this fraction was sufficiently large, the CMB and the mass power 
spectrum could be reproduced for the simple Chaplygin gas. To answer this question quan- 
titatively, we follow the Press-Schechter procedure [58] as in [38]. Assuming 5R(ad ec ) is given 
by a Gaussian random field with dispersion o~(R), and including the notorious factor of 2, to 
account for the cloud in cloud problem, the condensate fraction at a scale R is given by 

1 ' Js c (R) V2^a(R) P V 2a\R)) \y/2 a(R) J 1 ' 

where 6 C (R) is the threshold shown in figure [2j In figure [2] we also exhibit the dispersion 

a\R) = / — exp(-k 2 R 2 )A 2 (k, a dec ), (72) 
J o k 

calculated using the Gaussian window function and the variance of the concordance model 
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Figure 3: Fraction of the tachyon gas in collapsed objects using 5 C (R) and cr(R) from figure 
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A 2 (k,a) 



const 




T 2 (k) 



k 



n s — 1 



(73) 
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In figure 3 we present F(R), calculated using fl7TT) - (!73j) with const=7.11xl0~ 9 , the spectral 
index n s =1.02, and the parameterization of Bardeen et al [59] for the transfer function T(k) 
with f2^=0.04. The parameters are fixed by fitting (1731) to the 2dFGRS power spectrum 
data [60J. Our result demonstrates that the collapse fraction is about 70% for n = 2 for a 
wide range of the comoving size R and peaks at about 45% for n — 1. 

Albeit encouraging, these preliminary results do not in themselves demonstrate that the 
tachyon with potential (8) constitutes a viable cosmology. Such a step requires the inclusion 
of baryons and comparison with the full cosmological data, much of which obtains at low 
redshift. What has been shown is that it is not valid in an adiabatic model to simply pursue 
linear perturbations to the original background : the system evolves nonlinearly into a mixed 
system of gravitational condensate and residual k-essence so that the "background" at low 
z is quite different from the initial one. Because of this one needs new computational tools 
for a meaningful confrontation with the data. 



15 



5 Summary and Conclusions 



The first key test for any proposed quartessence model should be: Does it actually yield 
nonlinear dark matter structure, as well as linear dark energy in the inhomogeneous almost 
FRW universe that we see. In this paper we have analyzed the nonlinear evolution of the 
tachyon-like k-essence with a very simple potential V = Vop 2n . We have demonstrated that 
a significant fraction of the fluid, in particular for n = 2, collapses into condensate objects 
that play the role of cold dark matter. No dimensionless fine tunings were required beyond 
the inevitable \/~GH -C 1. 

Moreover, these results were obtained in a relativistic framework for nonlinear evolution 
that is as simple as the spherical dust model but includes the key effects of the acoustic 
horizon. Although it could be subject to possible improvements (e.g. a variable Gaussian 
width [38], and it is lacking the sophistication of the exact spherical model, [61], it does allow 
us to make the sort of quantitative assessments that have been missing [H2]- It is directly 
formulated in a convenient coordinate gauge, does not involve any hidden assumptions, and 
it easily deals with multi-component systems. 

The tachyon k-essence unification remains to be tested against large-scale structure and 
CMB observations. However, we maintain, contrary to the opinion advocated in [21], that 
the sound speed problem may be alleviated in unified models no more unnatural than the 
ACDM model. Indeed, an encouraging feature of the positive power-law potential is that 
it provides for acceleration as a periodic transient phenomenon [63J which obviates the de 
Sitter horizon problem [41 J , and, we speculate, could even be linked to inflation. 



A Adiabatic speed of sound 

The standard definition of the adiabatic speed of sound is 

2 _ 

Cs dp ' 

' s/n 

where the differentiation is taken at constant s/n, i.e. for an isentropic process. Here 
s = S/V is the entropy density and n = N/V the particle number density associated with 
the particle number N. We use the terminology and notation of Landau and Lifshitz [Bl] 
(see also [65])). For a general k-essence, with C = C(p,X), equation (f?4"l) may be written as 



(74) 



dp 
dp 



{dp/dX)dX + (dp/d<p)dip 



s/n 



(dp/dX)dX + (dp/d(p)dcp 



s/n 



where the differentials dX and dip are subject to the constraint d(s/n) 
that this constraint implies dip = 0. 

We start from the standard thermodynamical relation 



(75) 



0. Next we show 



d( p V) = TdS-pdV, 

where the volume is, up to a constant factor, given by V 
written in the form 

dh = Td I - ) + -dp , 



(76) 

1/n. Equation ( 1761) may then be 

(77) 
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where 

h = P -±l (78) 
n 

is the enthalpy per particle. For an isentropic relativistic flow one can define a flow potential 
if such that 

hu^ = (p t(t . (79) 

Comparing this with (TTTT) we find 

h = Vx. (80) 
This together with f|T2|) and (|T3|) yields in turn 

n = 2VxC X - (81) 

This expression for the particle number density is derived for an isentropic process. In a 
purely kinetic k-essence with C = C(X), equation flHTl) follows from the field equation for ip 

(£x<rV, M ) ;i , = 0, (82) 

which implies conservation of the current 

= 2£ x <p, fl = nu^ . (83) 

The particle number density n in this expression coincides with flHTl) . However, in a general 
k-essence, with C = C((p,X), the field equation (jTSj) does not imply that the current (1831 
is conserved. Nevertheless, equation (I8~T|) is still a valid expression for a conserved particle 
number density when the condition d(s/n) = is imposed. 
From (1771) with d(s/n) = and using (IHUj) we obtain 

ii 

dp = ndh= —=dX. (84) 
2 v X 



Comparing this with the general expression for the total differential of p 

d ^ = W dx + T^ (85) 

we must have dip = and 

d P n r fQR\ 

m=57! m 

as it should be. Hence, we conclude that an isentropic process implies dip = and equation 
r5] l yields 

2 #P 



° s dp 



P /dX 



dp/dX 
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B Braneworld Connection 



It is useful to view the tachyon condensate from the braneworld perspective. Consider a 3+1 
brane moving in a 4+1 bulk spacetime with metric 

dsl = g {5)MN dX M dX N = f(y) 2 g^(x)dx»dx v - dy 2 (88) 

where we generalize [9l [Hll H3] to allow for a warping of the constant y slices through f(y). 
The points on the brane are parameterized by X fl (x fl ), and G^ v = g^MN^^X^ is the 
induced metric. Taking the Gaussian normal parameterization X M = (x M ,Y(x M )), we have 

G, u = f(Y) 2 g^(x)-Y^ u . (89) 

The Dirac-Born-Infeld action for the brane is 



Sbrane = —O J d*XyJ- det G ^ = -O J d^Xy^g /{Yf 

and with the redefinitions 



f 2 (Y) 



(90) 



jfej = <P# , °f(Y) 4 = (91) 



we obtain 



^brane = -V(tp) yjl - g^^u • (92) 

For an unwarped bulk we obtain / = 1 and V = a = y/~A, i.e. the simple Chaplygin gas. 
In general, V(cp) identifies £brane with the tachyon model. If the brane also couples to, e.g., 
bulk form fields, there are additional terms that are functions of <p only. Thus every tachyon 
condensate model can be interpreted as a 3+1 brane moving in a 4+1 bulk. Note that this is 
true only for fl9"2l or (j521) but not for flBTl) . The prescription (1591 does not take into account 
the distortion of the bulk metric when the brane is not flat. This, however, can be accounted 
for using the methods of [66] . 

Given V(ip) the warp factor can be reconstructed via 

Y-Y = J f{Y{^))dv = a" 1 / 4 J V(ip)^*Ap. (93) 
For example, for the power-law potential 

V((p) = V oV 2n , (94) 

we obtain 

'M-rO • l -— 2 k) ■ (95) 

Using (l67j) and writing \/ Go = e/l, one finds e ~ (lHo) 1+n which is the only fine tuning for 
and a small-scale extra dimension. 
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